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ABSTRACT 

o : 

We discuss the observational properties of a remarkably faint triply-imaged 
galaxy revealed in a deep ^'-band Advanced Camera for Surveys observation of 
D . the lensing cluster Abell 2218 (z =0.175). A well-constrained mass model for the 

cluster, which incorporates the outcome of recent Keck spectroscopic campaigns, 
suggests that the triple system arises via a high redshift (z > 6) source viewed 
at high magnification (~ x25). Optical and infrared photometry from Hub- 
ble Space Telescope and the Keck Observatory confirms the lensing hypothesis 
0\ , and suggests a significant discontinuity occurs in the spectral energy distribu- 

^ , tion within the wavelength interval 9250-9850A. If this break is associated with 

Gunn-Peterson absorption from neutral hydrogen, a redshift of 6.6 < z <7.1 is 
inferred. Deep Keck spectroscopy conducted using both optical and infrared spec- 
trographs fails to reveal any prominent emission lines in this region. However, 
an infrared stellar continuum is detected whose decline below 9800A suggests a 
spectroscopic redshift towards the upper end of the range constrained photomet- 
rically, i.e. z ~7. Regardless of the precise redshift, the source is remarkably 
c3 ! compact (;$1 h^kpc) and faint (zfssolp =28.0) yet is undergoing vigorous star 

j> \ formation at a rate ~2.6 M yr _1 . An intriguing property is the steep slope of 

the ultraviolet continuum implied by the photometry which may suggest that the 
source is representative of an early population of galaxies responsible for cosmic 
reionization. Independent verification of these results is highly desirable but our 
attempts highlight the difficulty of studying such sources with present facilities 
and the challenges faced in pushing back the frontiers of the observable universe 
beyond z ~6.5. 
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Subject headings: cosmology: observations, galaxies: formation, galaxies: evolu- 
tion, gravitational lensing 

1. Introduction 

Reionization was a landmark event which imprinted a signature over the scale of the 
entire universe. After decades of lower limits on the redshift at which it occurred, recent 
observations of QSOs discovered by the Sloan Digital Sky Survey (Becker et al. 2001; Djor- 
govski et al. 2001, Fan et al. 2002) suggest that reionization was just finishing at z ~ 6 — 6.5. 
The discovery of z ~ 6.5 galaxies with strong Lyman a emission (Hu et al. 2002, Kodaira 
et al. 2003) is illustrative of possible sources which may be responsible. Analysis of recent 
temperature and polarization fluctuation data from the WMAP satellite suggests an optical 
depth for Thompson scattering of r = 0.17±0.04, implying that reionization began at higher 
redshift, perhaps as early as z ~ 15 — 20 (Kogut et al. 2003, Spergel et al. 2003). 

The discovery of star-forming galaxies at z ~ 6.5 is an important step toward under- 
standing the nature of the sources responsible for the end of cosmic reionization. However, 
to explore the earlier stages implied by the WMAP results, it is necessary to push the search 
for star forming systems to higher redshifts. With current facilities this is technically very 
challenging. It is perhaps salutary to note that the current redshift frontier (z = 6.5, corre- 
sponding to an observed Lyman a wavelength of 9200A) is coincident with the wavelength 
at which optical CCD detectors fall significantly in their quantum efficiency Ground-based 
infrared spectroscopy, necessary for exploring sources at higher redshift, is especially difficult 
at faint limits. 

Color-based searches for z >6 sources with the Advanced Camera for Surveys (ACS) 
on-board the Hubble Space Telescope (HST) is now recognized as a valuable way of locating 
z > 6 sources. Promising results have been obtained by utilizing the long wavelength F850LP 
(V-band) filter in conjunction with deep infrared imaging on HST 'or with large ground-based 
telescopes (Bouwens et al. 2003, Yan et al. 2003, Stanway et al. 2003, Dickinson et al. 2003). 
In view of significant contamination of red 'drop-out' samples by cool stars (e.g. Stanway 
et al. 2003), the primary challenge lies in spectroscopically verifying faint candidates (e.g. 
Dickinson et al. 2000). Neither optical nor infrared spectrographs on the current generation 
of ground-based telescopes may have the sensitivity to give convincing results unless strong 
emission lines are present. Although most of the distant sources found beyond z ~5 have 
been identified via strong Lyman a emission (Stern & Spinrad 1999, Spinrad 2003), some 
sources should have weak or no Lyman a emission (e.g. Spinrad et al. 1998). Indeed, a 
significant fraction (~ 75%) of the most intensely star-forming galaxies located by color 
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selection techniques at z ~ 3, reveal Lyman a only in absorption (Shapley et al. 2003). 

Gravitational magnification by foreground clusters of galaxies, whose mass distributions 
are tightly constrained by arcs and multiple images of known redshift, has already provided 
new information on the abundance of high redshift objects (Kneib et al 1996, Santos et al. 
2003). Particularly high magnifications (~ x 10-50) are expected in the critical regions which 
can be located precisely in well-understood clusters for sources occupying specific redshift 
ranges (Ellis et al. 2001). Although the volumes probed in this way are far smaller than 
those addressed in panoramic narrow band surveys (Hu et al. 1998, Malhotra et al. 2001, 
Hu et al. 2003) or the color-based surveys cited earlier, if the surface density of background 
sources is sufficient, lensing may provide the necessary boost for securing the first glimpse 
of young cosmic sources beyond z ~6.5 (Santos et al. 2003). 

In the course of studying the detailed rest-frame properties of the image pair of a lensed 
z = 5.576 galaxy in the cluster Abell 2218 (Ellis et al 2001), we have discovered a new faint 
pair of images in a deep ^'-band ACS observation of this cluster. Based on the geometrical 
configuration of this pair and its photometric properties, we concluded that the images most 
likely arise via strong magnification of a distant z >6 source. The implied high redshift of 
this source led us to explore its properties in more detail. 

A plan of the paper follows. We present the photometric observations in Section 2. 
Section 3 discusses redshift constraints determined independently from the lensing model of 
Abell 2218 and the spectral energy distribution based on broad-band photometry. Section 4 
summarizes our attempts to detect Lyman a emission spectroscopically and discusses the im- 
plications of a continuum discontinuity seen in the infrared spectrum. We discuss the source 
properties and implications further in Section 5. Throughout we assume a cosmological 
model with fi M = 0.3, fi A = 0.7 and H =70 km s" 1 Mpc" 1 . 

2. Photometric Observations 

The source in question was originally discovered as a pair of images (a and b on Figure 1) 
with reflection symmetry and a separation of 7" in a deep (5-orbit, 11.31 ksec) HST/ ACS- 
F850LP (z'-band) observation of the Abell 2218 cluster of galaxies. The observation was 
conducted as part of GO program 9452 (PI: Ellis) to characterize the stellar continuum in 
the lensed pair at £=5.576 discussed earlier by Ellis et al. (2001). 

The new ACS images were reduced using standard IRAF 4 and STSDAS routines. The 
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source was found by blinking the reduced F850LP image with two archival images taken at 
shorter wavelength with the F814W and F606W filters of the Wide Field Planetary Camera 2 
( WFPC2) (SM-3a ERO program 8500, PI: Fruchter; exposure time F606W: 10 ksec, F814W: 
12 ksec). A prominent discontinuity in the brightness of the pair can be seen viewing, in 
sequence, the F850LP, F814W and F606W frames (Figure 1). Together with the geometrical 
arrangement and symmetry of the pair in the context of our mass model for Abell 2218, this 
suggested the source is at high redshift and gravitationally magnified by the cluster (see §3). 

The two images have similar appearance in the HST/NICMOS F160W image. This H- 
band data was also obtained under the GO program 9452 and was acquired in a 4-orbit CVZ 
observation totalizing 22.2 ksec of exposure time. Similarly, we also identified the lensed 
pair in a deep Keck NIRC J-band image (5.64 ksec exposure) taken in 0.65"seeing on July 
22 and 23, 2002 (Blain and Reddy private communication). 

A clear mirror symmetry is seen in the HST data; both images contain a bright core, a 
second fainter knot and extended emission of lower surface brightness. The orientation of the 
pair also closely matches the predicted shear direction (Figure 1). The lensing hypothesis 
is further verified by comparing the colors of the two images a and b using the available 
photometric data summarized in Table l 5 . Colors were computed using a fixed elliptical 
aperture, and are identical for both images within the uncertainties. (For the z^lp — J 
color, the ACS z'-band image was convolved by a gaussian to match the seeing of the Keck 
NIRC J image). 



Association of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 

5 All photometric quantities are based on the Vega system 
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Figure 1: WFPC2-FQ06W, WFPC2-F8UW, ACS-F850LP and NICMOS-F160W images 
of Abell 2218 of the new faint pair in the lensing cluster Abell 2218 (2=0.175). The signals 
redward of the WFPC2-F814W observation suggests a marked break occurs in the continuum 
signal at around 9600A . Red lines correspond to the predicted location of the critical lines 
at z s =5,6.5 and 7 (from bottom to top, the latter two being almost coincident). The scale 
bar at the top left of each image represents 2". The predicted shear direction (thin blue 
lines) closely matches the orientation of the lensed images. 
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see jpeg file 

Figure 2: (Top) Location of the image pair a, b and the third image c in a pseudo-color 
image made from the WFPC2-F606W, F814W and ACS-F850LP images. The red curves 
refers to the critical lines of infinite magnification for sources placed at z=5.57Q and z=7.0 in 
the context of Kneib et al's (1996) mass model revised to include the ,2=5.576 pair (shown as 
unlabelled circles at the top of the figure) discussed by Ellis et al. (2001) and a triply-imaged 
SCUBA source at z=2.515 (Kneib et al. 2004a). (Bottom) Pseudo-color representation of 
the three images demonstrating their association with a single lensed source. 
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3. Redshift Constraints 

3.1. Gravitational Lensing 

In the context of the tightly-constrained mass model of Abell 2218 (Kneib et al. 1996), 
updated to include the properties of the ,2=5.576 pair identified by Ellis et al. (2001) and 
the more recent confirmation of a triply-imaged sub-millimeter selected source at z=2.515 
(Kneib et al. 2004a), the symmetry expected for a lensed pair around the critical line implies 
a source redshift z s >6 (see curves in Figure 1). The absolute location of the z=6 critical line 
is particularly well understood in this region from the measured symmetry of the adjacent 
2=5.576 pair around its critical line shown in Figure 2. However, the lensing configuration 
for the new source provides only a fairly weak constraint on the precise redshift beyond this 
lower limit since the location of the critical line does not change significantly beyond z ~6 
(Figure 1). 

Our mass model for Abell 2218 requires there to be a third image of the source, which we 
successfully located in the z'-band image at the expected position (image c in Figure 2) and 
with the expected flux (Table 1). Although our photometric coverage of this third, fainter, 
image is not as complete as that for the primary pair (because of the smaller field of the 
infrared cameras used), importantly the HST photometry confirms the same discontinuity 
in flux seen between F850LP and F606W (Figure 2, Table 1). The improved mass model 
suggests that all 3 images represent manifestations of a single source at z >6 magnified by 
a factor of ~25 (in the case of images a and b) . The intrinsic (unlensed) source brightness is 
ZS50LP = 28.0 ± 0.1, H imw = 26.5 ± 0.1. 

3.2. Spectral Energy distribution 

Figure 3 summarizes the available broad-band photometry for the brightest of the three 
images (a in Figure 1). We include, for completeness, the Keck NIRC J measurement 
although, as a ground-based measurement it is more adversely affected by crowding and 
background issues related to the adjacent luminous cluster members. Its low significance 
provides little more than evidence for a detection. Accordingly, we do not use the J band 
data in any of the subsequent analysis. 

A significant discontinuity in flux is apparent in the wavelength interval AA ~9200A- 
1/zm. The overlap in sensitivity between the WFPC-2 filter F814W and the ACS filter 
F850LP is of particular diagnostic use. If it is assumed the z > 6 source has a UV continuum 
which rises to shorter wavelengths with a discontinuity produced by a Gunn-Peterson trough 
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at around A rest =1216A, the ratio of the F814W and F850LP fluxes can be used to estimate 
the redshift depending on the slope of the UV continuum. 

Assuming the spectral energy distribution is given by a simple relation: /(A) oc \~ a for 
A/(l + z) >1216A and /(A) = otherwise, the accurate HST photometry (F606W, F814W, 
F850LP, F160W) implies the discontinuity occurs in the wavelength interval 9250-9850A, 
corresponding to 6.6< z <7.1 (see dashed lines in Figure 3) 6 . The photometric redshift 
constraint is particularly firm at the lower end. Below z ~6.6 it is hard to justify the observed 
F814W and F850LP flux ratio regardless of the form of the spectral energy distribution. In 
summary, therefore, the available HST photometry suggests the lensed source lies beyond 
z ~6.6. 



4. Optical and Infrared Spectroscopy 

Given the possibility that the lensed source lies beyond a redshift z ~6.6 with a Gunn- 
Peterson discontinuity in the wavelength range 9250-9850A we next tried to detect Lyman a 
emission in this region using both the Near Infrared Spectrograph (NIRSPEC, McLean et al. 
2001) on Keck II and the Low Resolution Imaging Spectrograph (LRIS, Oke et al. 1998) on 
Keck I. If around half of the flux in the F850LP band arises via a Lyman a emission line, in 
a manner analogous to the z =5.576 source (Ellis et al. 2001), we considered that the LRIS 
campaign should be successful in securing the redshift, particularly since the lower part of 
the region to explore is, by good fortune, one of the "windows" of low OH sky emission used 
by narrow band imagers to locate Lyman a emitters (e.g. Kodaira et al. 2003). 

We used LRIS on two runs, May 31 - June 1 2003 and June 30 - July 1 2003 (see 
Kneib et al. 2004b for a complete description of these data). During the first observing 
run we observed the brighter pair of the triple system using the 600 line grating blazed 
at lyum for a total of 13.8 ksec. The first night offered relatively poor conditions (cirrus, 
high humidity) but 9.0 ksec was secured on the second night in slightly better conditions. 
For these observations the wavelength coverage ranged from 6930 to 9500A. No signal was 
detected from either image. 

During the second observing run, we used the same configuration and integrated for 
9.2 ksec in relatively good conditions with ~ 0.8" seeing. The wavelength coverage was 
extended slightly to the red to reach 9600A. In the following, only the second observing run 



6 Although the weak NIRC J detection suggests a shallower UV slope, we give this discrepancy low weight 
in view of the superior quality of the HST data 
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Figure 3: (Top) Relative efficiencies of the filter+instrument used in our photometric obser- 
vations. From left to right: WFPC2-F8UW, ACS-F850LP and NICMOS-F160W. (Bottom) 
Spectral energy distribution of image a uncorrected for lensing magnification. Photometric 
points are indicated by blue crosses, and the 3a point source detection limits are indicated 
by a horizontal dashed red lines below each data point. The non-detection of the contin- 
uum in the LRIS 9000-9300A window is indicated by the black arrow. Red dashed lines 
correspond to power law spectral energy distributions (with / oc \~ a for A/(l + z) > 12 16 A 
and / = otherwise) with slope indices (from bottom to top at 1.2 10 4 A) of a=3,4 and 5. 
The available data are consistent with a neutral hydrogen break in the interval 9250-9850A 
corresponding to 6.6< z <7.1. 
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will be utilized as the data quality is much better. Flux calibration was conducted using 
Feige 67, Feige 110 and Wolf 1346 spectrophotometric standards observed in twilight. The 
spectroscopic observations were reduced using IRAF. Despite the improved conditions, no 
continuum or emission lines were detected to 9600A from either image. The absence of any 
LRIS detection in the clean wavelength range 9000-9300 A places an additional constraint 
on the lower redshift limit. In Figure 3 we illustrate this via an upper limit on the absolute 
flux (black symbol) which further suggests the source lies beyond z ~6.6. 

The NIRSPEC (McLean et al. 1998) observations were obtained with the Keck II tele- 
scope on the nights of May 10 and 11 2003. The spectrograph was used in low- dispersion 
mode using a filter with transmission from 9500 to 11200A. A slit width of 0.76" was used 
to acquire the pair giving a spectral resolution of i?=1100. We obtained 37 exposures of 
900 sec yielding a total integration time of 33 ksec. We dithered along the slit between two 
different positions, verifying our pointing on a nearby reference star using the slit-viewing 
guide camera at each dither. These observations were also reduced using IRAF. Individual 
spectra were registered using offsets determined from simultaneous images obtained with 
the slit-viewing guide-camera. Our optimally combined spectrum of each image comprised 
26 individual exposures totaling an effective exposure time of 23 ksec. All of the combined 
exposures were obtained in photometric conditions or through thin clouds. The data were 
flux calibrated using observations of Feige 110 (Massey & Gronwall 1990) which provide an 
effective calibration from 9700-10200A . Tests with two other standard stars (Feige 34 and 
Wolf 1346) indicate the relative flux calibration should be much more reliable than the ab- 
solute calibration; we have probably suffered from slit losses and possible absorption during 
the night due to faint cirrus. No emission lines were seen in the extracted NIRSPEC spectra 
but a faint stellar continuum was detected for both images a and b; the signal is slightly 
stronger for a. 

Although observational selection plays an important role, as most successfully identified 
z > 5 galaxies display intense Lyman a emission (Spinrad 2003), it is interesting to consider 
the maximum possible observed equivalent width, W max , implied by our non-detections in 
the LRIS and NIRSPEC data. Using our spectrophotometric calibration, assuming a line 
width of 5A and a conservative 5<r detection limit, Figure 4 shows W max as a function of 
wavelength for our data. Although OH airglow emission precludes detection in a few small 
wavelength regions, for 60% of the important 9000A-9500A interval an observed equivalent 
width larger than 120 A (a value much less than for most high redshift star forming sources, 
Hu et al 2003) would have been detected. From 9000-9300A (6.4< z <6.65) and longward 
of 9550A {z >6.85) the constraint is tighter. Because of the much deeper exposures with 
NIRSPEC and the stronger continuum redward of 9800A, a more stringent 5(7 upper limit 
of W max <60A is derived for the 9550-1 HOOA window effectively ruling out any reasonable 
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level of emission in the redshift range 6.85< z <8.2. Thus it seems reasonable to deduce 
that, if the redshift is not in the 6.65< z <6.85 interval, Lyman a emission is either very 
weak or absent. It will be important in future, deeper, observations to rule out an emission 
line hiding in the OH forest at ~9500A corresponding to a source at z ~6.8. 

The absence of Lyman a emission in a distant source may seem surprising. However, 
there are many examples of luminous Lyman break galaxies at lower redshift with weak or 
no emission (Shapley et al. 2003). Indeed, those authors claim only 20-25% of star-forming 
examples show Lyman a emission sufficiently prominent to be classified as narrow band 
excess objects. 

The most important outcome of the long NIRSPEC exposure is the detection of a very 
faint continuum redward of 9800 A (Figure 5). A continuum signal is seen in both images. 
As it is stronger in image a, we will use this spectrum for the following analysis. Given the 
absence of Lyman a emission and the inference of a photometric break in the wavelength 
range 9250-9850A the key question is therefore the lowest wavelength at which a continuum 
signal can be seen in the NIRSPEC data. 

Inspection of the data reveals a drop in the continuum flux at 9800A shortward of 
which there is no reliably detected signal despite a robust photometric calibration (see upper 
panels of Figure 5). The OH spectrum is fairly clean in this region although there is some 
atmospheric absorption which could affect the calibration at 9500-9700A . Although the 
feature itself is only marginal, the absence of flux below 9800A seems significant when one 
considers the photometric data summarised in Figure 3. Specifically, if the UV continuum 
extended down to 9250A , as would be the case if the source were at ^=6.6, it is difficult 
to understand why the stellar continuum is not detected to shorter wavelengths (Figure 5). 
Although this weak feature is the only indicator in our exhaustive attempts to measure a 
spectroscopic redshift, if it is indeed the cause of the Gunn-Peterson edge inferred from the 
photometric data, a redshift of z ~7.05 is implied. In this case, the photometric redshift 
analysis discussed in S§3.2 would indicate a steep UV continuum slope of a=5. 
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Lyman— a Redshift 




Wavelength (Angstrom) 
Figure 4 Observed equivalent width detection limit (5a), W max , for Lyman a emission 
in the LRIS and NIRSPEC data assuming a line width of 5A . An emission line stronger 
than H / =120A (indicated by the horizontal line), corresponding to an integrated line flux 
of 1.6xl0 -18 ergs s -1 cm -2 A -1 , would have been seen for 60% of the wavelength range 
9000A-9500A. Constraints in the range 9000-9300A and beyond 9550A are tighter. The 
longer NIRSPEC integrations provide a maximum equivalent width below W=60A for the 
entire 9550-lllOOA range, corresponding to an integrated line flux of 7x10 ergs s 
cm~ 2 A _1 . 
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Figure 5 NIRSPEC spectrum for image a in the wavelength region 9500A-l.l/xm. (Top) 
Relative sensitivity of NIRSPEC in the configuration used derived from 3 independent flux 
calibration stars. Adequate sensitivity is available over the entire range. (Middle) observed 
spectrum of one of these standards, Feige 110 illustrating a priori the absence of any strong 
atmospheric features in the reduction. The red hashed rectangle indicates the region where 
variable atmospheric absorption is expected (Massey & Gronwall 1990). (Bottom) Flux cal- 
ibrated binned spectrum for image a. Error bars include the effect of wavelength-dependent 
OH emission. Curves represent UV SEDs consistent with the photometry (Figure 3) span- 
ning the range 6.6< z <7.1. The absence of significant flux below 9800 A suggests a redshift 
at the upper end of this range, i.e. z ~7. 
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5. A z~7 Source 

Clearly identifying the redshift of this source has been extremely challenging and further 
studies would be highly desirable to confirm our conclusions. We have introduced three 
independent arguments which, taken in combination, justify why the newly-located multiple 
images arise via a single lensed source at a redshift z ~7. First, the geometrical configuration 
of images a, b and c with respect to the critical line in this well-modeled cluster indicates 
z >6. We argued that the location of the critical line in the vicinity of the pair is tightly 
constrained from the successful identification of the earlier source at £=5.576 (Ellis et al 
2001). Secondly, there is a strong break indicated in the spectral energy distribution (SED) 
delineated by the broad-band photometry extending from the F606W to F160W (H) filters, 
arguing for a redshift 6.6 < z < 7.1. The lower redshift limit z >6.6 is particularly firm and 
supported by the absence of any LRIS continuum signal in the wavelength range 9000-9300 
A. Finally the absence of any detectable signal below 9800A seen in the NIRSPEC continuum 
argues the redshift lies in the upper end of the photometric range, z ~7, as otherwise the 
rising UV continuum would have been detected. 

Even if one disregards the weak feature in the NIRSPEC spectrum, the 6.6< z < 7.1 
source is of considerable interest since, as we discussed in §3-4, it is difficult to reconcile a 
spectral break at A >9250A with a UV continuum slope of a^3, assuming the UV continuum 
is described by a power law above the Lyman a discontinuity (Figure 3). Starburst models 
assuming Population II metallicity and a Salpeter stellar initial mass function (IMF) typically 
produce slopes of a^2 (Leitherer et al. 1999). Clearly the inferred UV continuum of the z ~ 7 
source rises faster to shorter wavelengths than for any of the models calibrated with local 
data. Models based on massive metal-free stars, such as those which may be expected close 
to the epoch of reionization (Abel, Bryan, & Norman 2000; Bromm, Coppi, & Larson 1999), 
do produce somewhat steeper UV slopes, a ~ 3 (Schaerer 2002). More precise photometry 
for this object, particularly in the 1-1.5/im range is needed to understand this issue. If a 
steep UV slope (a > 3) is indeed confirmed, the source may represent a promising candidate 
for a Population III starburst. 

The star-formation rate (SFR) of the source can be estimated from the stellar UV 
continuum luminosity. The unlensed F160W magnitude, H ieow = 26.5 ± 0.1, translates 
into a specific flux of 2.35 x 10" 31 erg s _1 cm~ 2 Hz _1 at 1.6/xm, assuming a smooth SED. The 
intrinsic specific luminosity at 2000A is then 1.85 x 10 28 ergs _1 Hz -1 implying a SFR of 
2.6 M yr _1 using Kennicutt's (1998) empirical calibration. We estimate a 15% error arising 
from uncertainties in the photometry and magnification. This derived SFR is significantly 
higher than that for the ~0.5 M & yr _1 computed for the lensed source at £=5.576 (Ellis et 
al. 2001) but within the range of typical Lyman-a emitters at z ~5.5-6.5. 
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Our estimated SFR may suffer from two systematic errors. If dust is present, it will dim 
the UV continuum, causing us to underestimate the value. However, assuming a familiar 
selective extinction curve, reddening would then imply the intrinsic UV continuum slope is 
even steeper than a ~3-5. A second uncertainty arises from our assumed stellar initial mass 
function (IMF) and metallicity via Kennicutt's local calibration. If the IMF is top-heavy or 
the metallicity is lower than in local starbursts, the SFR will likely have been over-estimated. 

Given the intense star formation, the apparent lack of Lyman a emission is puzzling if 
dust extinction is not important. Weak or absent Lyman a emission seen in the younger (<1 
Gyr) Lyman break galaxies at z ~3 has been interpreted via the presence of dust shrouds 
which are eventually disrupted via feedback processes as the stellar population matures 
(Shapley et al. 2003). As our source is likely younger than 500 Myr, a very specific dust/gas 
geometry would be needed to strongly extinct the Lyman a photons without reddening the 
UV continuum. This problem may be exacerbated if the IMF is top-heavy or the metallicity 
low, since the ionizing photon production rate will be higher making any emission line yet 
more prominent (c.f. Malhotra & Rhoads 2002). 

An alternative explanation for the absence of Lyman a emission may be incomplete 
reionization at z ~7. Neutral hydrogen can scatter Lyman a photons after they have escaped 
the ISM of the emitting galaxy (e.g., Miralda-Escude & Rees 1998). If the source were 
embedded in a neutral zone, even a strong Lyman a emission line emitted from the galaxy 
could be quenched. However, in such cases, simulations suggest that emission can still be 
observed even from sources embedded in fully neutral zones, depending on details of the 
many relevant parameters of the source (Haiman 2002, Santos 2003). 

Since the universe is only 750 Myr old at z — 7, the total stellar mass is unlikely to 
exceed 1 — 2 x 10 9 M . Halos of mass > 1.2 x 10 10 M Q at z ~7 are quite consistent with 
conventional structure formation models (Barkana & Loeb 2002) and sufficient to supply 
baryons to sustain the observed star formation rate. The instantaneous baryonic accretion 
rate for such a halo at z = 7 is, on average, 4.4 M yr _1 (e.g. Lacey & Cole 1993), so even if an 
earlier starburst consumed most of the baryons, they can be rapidly replenished. Although 
we have no constraints on the actual star-formation history, these arguments emphasize that 
we are not necessarily viewing the system at a special time in its evolution. 

The strong magnification (~ x25) of the two brighter images gives us our first glimpse 
into the morphological structure of a very distant source on sub-kpc scales. Both images 
a and b have a bright core, a second fainter knot, and extended emission of lower surface 
brightness (Figure 1). There is no noticeable color gradient. Our mass model for Abell 2218 
implies magnifications of ~ xl5 along the major axis of images a and b, and ~ xl.7 along 
the minor axis. The 0.15" width of the F850LP PSF translates into a physical resolution of 
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~470 pc along the minor axis and ~50 pc along the major axis. 

We thus infer that the source (observed with 3.6" x < 0.15") has a maximum physical 
size of 1.2 kpc by <500 pc. The maximal associated area of 0.6 kpc 2 indicates a star- formation 
surface density in excess of 4.3 M yr" 1 kpc" 2 . The bright knot is only 100 pc across implying 
a star-formation surface density in the range ~50-250 M yr _1 kpc _2 depending on the exact 
geometry, comparable to the most intense starburst activity observed locally (Kennicutt 
1998). 

If this source is typical of those which reionized the Universe in a narrow time interval 
of A z—1 around z ~7, we can estimate the expected surface density from the arguments 
presented by Stiavelli et al (2003). Depending on the source temperature, Lyman continuum 
escape fraction and dumpiness of the IGM, we would deduce surface densities n ~0.3-5 
arcmin -2 are necessary. In as much as it is possible to estimate the actual surface density 
from one source and the limited area examined by looking through only one cluster lens to 
find it, we find a number density of n ~ 1 ± 0.5 arcmin -2 . 

In conclusion, we present evidence that we have found a highly magnified source which 
lies beyond z ~ 6.6, possibly at z ~7. Even in advance of the infrared capabilities of JWST, 
further observation of this source will be important in determining a more secure redshift, 
and improving the constraints on the slope of the UV continuum. Further spectroscopy in 
the 9200A-l/xm region would be valuable to probe any Lyman a emission at z ~6.8 (Figure 
4) and to confirm (or otherwise) the significance of the weak continuum drop seen in the 
NIRSPEC data at 9800A . The location of the Gunn-Peterson trough might also be verified 
more precisely via narrow band imaging through gaps in OH forest and with a deep ACS 
grism spectroscopy. 

Notwithstanding the need for further work and regardless of its precise redshift in the 
constrained window 6.6< z <7.1, our source appears to be a star-forming galaxy with 
intriguing spectral properties, possibly representative of a new population responsible for 
ending cosmic reionization. While our discovery highlights many of the challenges facing 
searches for those z >6.5 galaxies responsible for reionization, it also demonstrates the 
ability of strong lensing by clusters of galaxies to locate and reveal the detailed properties 
of high redshift sources. The lensed galaxy presented in this paper, if observed unlensed 
at zab ~ 28.5 would lie at the detection limit of the upcoming UDF observation 7 and no 
spectroscopic follow-up would have been possible. 



7 http: / / www.stsci.edu/hst / udf/parameters 
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Table 1. Observed Photometry for the Triple System 





a 


b 


c 


a J2000 


16:35:54.73 


16:35:54.40 


16:35:48.92 


^J2000 


66:12:39.00 


66:12:32.80 


66:12:02.45 




not detected 


not detected 


not detected 


huw 


26.5±0.2 


26.4±0.2 


not detected 


^850LP 


24.38±0.05 


24.54±0.05 


25.9±0.1 


J 


23.4±0.3 


23.5±0.3 






22.96±0.06 


23.01±0.07 




V606W-Z850LP 


>3.6 


>3.5 


>2.1 


hl4W~ z 850LP 


2.1±0.2 


1.9±0.2 


>1.8 


Z850LP-J 


1.0±0.3 


1.0±0.4 




Z850LP-H\qow 


1.42±0.1 


1.53±0.1 




Magnification 


25±3 


25±3 


5.3±0.5 



3a detection limits for a point source: 

^606W — 28.0, Isuw — 27.2, ZgsoLP — 26.7, J = 24.4, H 1&0 w = 25.8. 
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